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ABSTRACT 

Correlation coefficients have been computed to 
investigate the parameters used to describe the spherical 
model of an eclipsing binary system. Regions in 
parameter hyperspace have been identified where strong 
correlations exist and, by implication, the solution 
determinacy is low. The results are presented in 
tabular form for a large number of system configurations. 

INTRODUCTION 

Russell and Merrill (1952) have pointed out that indeterminacy 
in the elements of eclipsing binaries can persist even though 
high quality observations are available and regardless of the 
formal method employed to solve for the elements. Such a 
situation exists for the case of shallow partial eclipses; large 
changes in the ratio of the radii result in only minor changes 
in the shapes of the eclipses. Another example is the 
low determinacy of the limb-darkening coefficients for 
partial eclipses. In Russell and Merrill’s discussion, 
generalization of the determinacy of a solution is based on the 
size of the angle formed at the intersection of the shape and depth 
relations. The shallower the intersection the 
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lower the determinacy. A number of troublesome situations 
are described by Russell and Merrill and these results will 
not be repeated here. Kopal (1959) has discussed this problem 
also and points out that the rate of convergence is related 
to the determinacy of the elements and that a solution depends 
upon the nature of the eclipses. Again, shallow partials 
produce low determinacy. 

It is difficult to express the determinacy in a closed 

analytic form and for this reason a numerical approach is 

used in the present work. Rather than considering the weight 

or determinacy of the parameters directly, it was found more 

convenient to employ the correlations among the parameters as 

indicators of the reliability of the elements. The primary 

objective of this work then is to identify those regions in 

parameter hyperspace where the correlations are most serious and, 

by implication, the determinacy the poorest. In practice one is 

required to use data from both eclipses for partials due to the 

ill-conditioning of the eclipse relations. To simulate this 

practical aspect, the correlation coefficients have been 

computed with model data from both eclipses. Preliminary 

results indicated that the maximum geometric eclipse depth, 

p is a suitable choice for the independent variable in 
o 

order to demonstrate the trends in the correlations in an 
efficient and meaningful way. The matrix of correlation 
coefficients is found for a range of geometric depths at grid 
points representing a wide range of eclipsing configurations. 
Only spherical stars with no radiative interactions are treated. 

Method 

A program for calculating the light and light loss during 
eclipses of spherical stars was written by the authors for the 
IBM 360 computer at GSFC. Although in its most general form 
eccentric orbits and non-linear limb-darkening can be treated, 
the eccentricity was set to zero and the non-linear terms for 
the darkening were suppressed in this work. In this initial 
treatment of the problem it was decided to express the model 
in terms of the most basic set of parameters. These are given 


below. 



g’ 


u 


e* 


R 


- radii of the larger and smaller stars, respectively, 

i 

in units of the semi-major axis. 

- inclination of the orbital plane relative to the 
plane of the sky. 

- linear limb-darkening coefficients of the respective 
stars . 

- ratio of surface brightnesses, the cooler to the 
hotter star. 

The light of the system is given by the following: 

/ - 1- i Jj I d a 

where the limits of the integral refer to the eclipsed portion 

of the eclipsed star, d^ is the elemental surface area, L is 

the sum of the light of the uneclipsed components, and 

I = I (1 - u + up) expresses the conventional linear limb-darkening 

relation. 

In an orthogonal system this equation can be integrated 

analytically over one dummy variable. Tests are made within 

the program to determine the limits for the first integration 

as a function of the geometry of the eclipse appropriate for the 

given mean anomaly. Then the resulting expression is evaluated 

numerically by Gaussian quadratures using double precision 

arithmetic throughout. A ten point Gaussian form is used. 

The accuracy is such to reproduce the tabulated a values given 

by Russell and Merrill. At each grid point, and maximum geometric 

depth the light is calculated for 40 equidistant phase values 

between minimum and external tangency for the ascending branch 

of both eclipses. Of course J = 1 outside of eclipses. 

The correlation coefficients, p^, °f the six parameter 

models are given in terms of the elements, E. of the moment 

~ 1/2 

matrix as P. . = E. ./(E. .E ) . In turn, the moment matrix is 

^ J J 1 J J 

defined as twice the inverse of the matrix of the second partial 

derivatives R. . where, in terms of the light, X, R. . ( 6,,...9_ 

ij ij J- o 

00. 00.. The 0.*s are the six parameters defined above. In order 
i J 

to compute the derivatives, a total of 73 distinct models must 



be computed for each grid point at each geometric depth; one 
to define the run of the light curve for the parameters given 
at the grid location and the other 72 for the light curve 
variations induced by the cross term perturbation of the parameters. 
Subsections of the subroutine STEPIT, written by J. Chandler (1965), 
were used to organize these computations and perform the matrix 
inversion. The results appear in Tables 1-5. Unless noted the 
computations are for r g = 0.25 and equal limb-darkenings , u = u g = 0.6, 
but Tables 4 and 5 show results for significant changes in these 
base parameters. The explanation of the tabulated material is 
as follows: 

Column 1 ratio of the radii, occultations have k > 1.0. 


2 

geometric depth, 
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A row of zero values means that no computations were made for 
that particular geometric depth. A negative sign means anti- 
correlation. 

Discussion 

No attempt will be made to discuss the form or significance 
of ■ all the correlation trends apparent in this material. or to 
generalize them. The intention of this work is to provide a 
framework within which the analysis of individual solutions can 
be assessed. High correlations imply poor convergence and 
certainly indicate anomalous error assignments for the solution 
parameters. For perspective, the commonly accepted rule is 
that p < J.5 implies no correlation while p>|o.8| implies high 
correlation. The generally conceded increase in the indeterminacy 
for shallow partial eclipses is evident and is quantified in 
this work. The high determinacy for total eclipses is verified 
except that a moderately strong correlation between u_ and R 

o 



appears to exist. As R approaches zero, the main parameters 

r . r . and i become strongly correlated. Troublesome 
g’ sj 

correlations also appear in partial occultations for internally 
tangent configurations. 

From a different point of view these results tend to reaffirm 
the assertion made by Kopal that a solution is not complete unless 
the errors of the elements are evaluated. As a future 
investigation, it would appear desirable to study selected 
combinations of parameters which would minimize the correlations 
or remove the skewness of the probability error ellipse and 
thereby improve the determinacy of the solution. 
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Table 1 


CORRELATIONS FOR A L 1 GET RATIO = 0.9 


Q m SC 

1 . 56 

C. 35 

-C, 43 

0. 63 

0. 08 

-0. 01 

-0.77 

-0.35 

o.sc 

1 . 52 

0.4-5 

— -0. 56 

0. 56 

0. 2 0 

0.10 

-0, 8 7 

-0 .29 

0.50 

1 .58 

Oc 37 

— Co 4 4 

0, 57 

0. 1 9 

0.05 

-0.89 

-0 . 35 

c • s c 

1 .45 

0.47 

-0.61 

0. 78 

0.4 5 

0. 1 1 

-0.90 

0.21 

0.50 

1 .40 

0.14 

-0. 35 

0. 83 

0.82 

0.49 

-0.94 

C .19 

o.sc 

1 . 36 

-0.68 

0,64 

0. 73 

0. 87 

0. 69 

-1.00 

-0 .06 

0.50 

1 . 32 

-0. 43 

0, 53 

0. 60 

0. 74 

0. 03 

-0. 99 

-0 .33 

0.5 0 

1 . 28 

C. 0 

0. 0 

0, 0 

0. 0 

0.0 

0.0 

0 .0 

0.P0 

1 .56 

-0. C7 

-0. 21 

0. 70 

-0. 01 

-0. 25 

-0.60 

-0 .37 

O.PC 

1 . 52 

0.12 

-0. 61 

0, 83 

0, 37 

-0. 21 

-0,6 3 

0.13 

0.80 

1 . 48 

0.17 

-0.54 

0. 74 

0. 62 

0.13 

-0.86 

0 .57 

O.SC 

1 . 44 

Co 25 

-0.63 

0, 62 

0, 84 

0. 44 

-0, 86 

0 .66 

0.6 0 

1.40 

0, 1 7 

-0. 7 3 

0. 90 

0. 92 

0. 54 

-0.73 

0.51 

O.SC 

1 . 36 

-0.43 

-C. 68 

C. 39 

0.92 

0.40 

-C.36 

-0.06 

o.eo 

1 . .12 

-0.65 

-0. 26 

0. 32 

o.e3 

0.03 

-0.50 

-0 .28 

0.50 

1.20 

-0.74 

C. 27 

0. 73 

0. 56 

-0.37 

-0.75 

-0 .34 

1.0C 

1.56 

-0.37 

-0.11 

0. 32 

-0.25 

-0.29 

-0.1 1 

-0.32 

1 .cc 

1 . 52 

-C. 61 

-0.30 

0. 34 

0. IS 

-0. 07 

-0.2 9 

0 .03 

1 . 00 

1.48 

-0.42 

-C.51 

0. 4 8 

0.50 

0.09 

-0.4 9 

0.39 

1.00 

1.44 

0.15 

- C. 75 

0. 74 

0.77 

0.12 

-C.74 

0 .72 

1 .CO 

1.40 

0.66 

-0.91 

C. 91 

C. 92 

0.03 

-0.9 1 

0 .90 

1 .00 

1.36 

C. 13 

-0. 77 

0. 76 

0.80 

0. 10 

-0.76 

0 .74 

1 .00 

1. 32 

-0.40 

-0.54 

0. 55 

0. 61 

0. 05 

-0.55 

0 .47 

l .oc 

1 • 28 

-0.62 

-0, 33 

0. 4 7 

0.44 

-0. 14 

-C.4S 

0.27 

1 o25 

1 . 56 

-Co 04 

— 0. 57 

0. 33 

-0. 39 

-0.43 

-0.29 

-O.C1 

1.25 

1 . 52 

0. 1 3 

-0. 62 

0, 38 

0. 1 1 

-0.23 

-C.6 1 

0 . 39 

1.25 

1 . 46 

0.16 

-0.85 

0. 5 3 

0.53 

-0. 16 

-0.54 

0 .53 

1.25 

1 . 44 

0.34 

-0.83 

0. 7C 

0. 73 

-0. 26 

-0.73 

0 .33 

1.25 

1.40 

C. 39 

-C.0O 

0. 65 

0.69 

-0. 26 

-0. 86 

0.93 

1.25 

1 . 36 

-C. 26 

-0.44 

0. 13 

0. 1 9 

0.04 

-0.75 

0 .89 

1.25 

1.22 

-0.62 

— 0. 45 

-0. 05 

-0. 07 

-0. 10 

-0.40 

0.75 

1.25 

1 . 28 

-C. 72 

-C. 67 

0. 1 3 

-0.23 

-0.47 

0. 1 2 

0 .4 2 

2.0C 

1 . 56 

0 > 28 

-0. 75 

C, 3 5 

-0. 38 

-0. 38 

-0.50 

0.10 

2. CO 

1 . 52 

0.55 

-0.86 

0. 5 3 

-0. 30 

-0. 52 

-0, 6 4 

0 .29 

2.0 0 

1.46 

C . 45 

-0, 57 

0, 51 

-0. 34 

-0.51 

-0.5 5 

0 .25 

2.00 

1 . 44 

0 c 5 1 

-C. 88 

0. 3 2 

0. 20 

-0.21 

-0.67 

0 .45 

2.CC 

1.40 

0.25 

-C, 91 

0, 2 7 

0.26 

-0.13 

-0. 5 2 

0 .84 

2.0 0 

1 . 36 

- 0 o 6 3 

-0. 9 3 

-0. 2 1 

-0. 02 

0.4i 

0.55 

0 ©67 

2 « C 0 

1 e 32 

- 0 4 3 

-0, 99 

-C, 47 

-0. 34 

-0. 02 

0.9 2 

0 o75 

2 . C C 

1 . 26 

0 3 C 

0,0 

c, c 

0. C 

0. 0 

C.O 

0 .0 

( 1 ) 

( 2) 

( 3 ) 

( 4 » 

( 5 ) 

(6 > 

(71 

(3 I 

(9) 


0.25 

0.25 0.03 -0.31 

- 0. 2 8 

- 0. 09 

-0. 21 

0, 99 

0. 37 

0.36 0.05 -0. 42 

- 0. 39 

- 0. C9 

-0. 22 

0, 99 

0.36 

0.37 0. 14 -0.42 

- 0. 40 

-0.07 

-0. 27 

0, 97 

0. 25 

0.1 2 -0. 40 -0.37 

- 0. 1 8 

0. 36 

-C. 10 

0, 89 

0. 14 

0.02 -0.29 -0. 32 

- 0.2 1 

0.72 

0.20 

0. 92 

-0.28 

-0.48 0.C4 0.24 

0.4 2 

0. 90 

0. 33 

0. 71 

-0.47 

— 0 .02 0.39 0.50 

- 0. 0 5 

0.90 

-0.34 

0. 09 

0.0 

0.0 0.0 0.0 

0.0 

0. 0 

0. 0 

0, 0 

0. 30 

0.36 -0. 23 -0.43 

- 0.26 

-0. 11 

-0. 45 

0. 92 

0.3 1 

0.16 -0.76 -0.57 

0.02 

0.37 

-0. 35 

0, 73 

0.5C 

0.08 -0.86 -0. 74 

- 0. 1 4 

0.65 

-0. 09 

0, 69 

C.53 

0.15 -0.92 -0. 88 

- 0. 3 5 

0. 87 

C. 24 

0, 67 

0.46 

0.10 -0.94 -0.92 

■ 0. 4 4 

0. 94 

0. 37 

0, 66 

-0. 1 3 

— C .16 - 0.86 -o.e«- 

- 0. 3 1 

0. 93 

C. 20 

0, 54 

-0. 30 

0.C5 -0.50 -0.55 

- 0. 1 9 

o. ee 

-0.20 

0, 32 

-0. 06 

0.46 C. 22 -0. 1 1 

- 0.4 6 

0.67 

-0. 54 

0, 24 

0. 19 

0.26 -0.27 -0.23 

-0.01 

-0.37 

-0.80 

0, 84 

C. 18 

0.05 -0.68 -0.65 

-0. 0 1 

0. 15 

-0.61 

0, 67 

0.33 

-0.10 -0.89 -0.86 

0.0 

C. 65 

-0. 38 

0, 45 

0.66 

-0.13 -0.97 -0.96 

0.0 

0.91 

-C. 18 

0, 23 

0.83 

-0.08 -0.99 -0.99 

0.0 

0.98 

-0. 08 

0. 1 1 

0.69 

-0 . 10 -0. 93 -0. 97 

0.0 

0.95 

-0. 13 

0, 16 

0. 39 

-0.04 -0.92 -0.90 

0. 3 

0.87 

-0. 22 

0. 29 

0.29 

0.15 -0.74 -0.72 

0. 0 

0. 76 

-0. 32 

0, 36 

0.64 

0.25 -0.48 -0.26 

o, 2 e 

- 0. 16 

-0. 91 

0, 52 

0. 70 

0.22 -C.59 -0.76 

0.0 

0. 32 

-C. 72 

0, 41 

0.68 

-0.12 -0. 79 -0.86 

0. 1 7 

0.63 

-0. 67 

0, 14 

0.80 

-0.44 -0.92 -0.94 

0.4 1 

0.88 

-0.66 

-0, 23 

0.9C 

-0.54 -0.96 -0.96 

0. 5 0 

0.95 

-0.64 

-0. 38 

0.85 

-0.31 -0.91 -0.91 

0.29 

0. 92 

-0.46 

-0. 07 

0.72 

0.04 -0.72 -0.62 

0. 1 9 

o. ei 

-C. 27 

0, 34 

0. 64 

0.42 -0.29 0.13 

0.47 

0.59 

-0. 25 

0. 62 

0.56 

-0.03 -0.43 0.02 

0. 4 1 

-0.20 

-0.99 

0. 30 

0.48 

-0.21 -0.58 0.03 

0.5 4 

-0.19 

-0.99 

0, 30 

0.49 

-0.15 -0.57 0.C9 

0.5 4 

- 0. 19 

- C. 98 

0, 35 

0.7 1 

-0.16 -0.45 -0.40 

0. 2 3 

0.22 

-C. 91 

0, 19 

0,76 

-0.56 -C.E1 -0.38 

0. 39 

0. 64 

-C. 84 

-0„ 13 

0.73 

-0.67 0.12 -0.04 

- 0,3 1 

0.91 

-Co 67 

-0, 30 

0. 6 1 

0 .0 1 0.49 0.6 1 

0. 1 3 

0.91 

-C. 02 

0, 39 

0.0 

0 .0 0.0 0. 0 

0. 0 

0.0 

0. 0 

0. 0 

(10 ) 

(11) (12) (13) 

( 14 ) 

( 15) 

< 16) 

(17) 


See text for explanation of tabular material 



Table 2 


CORRELATIONS for a light ratio = 0.6 


0.50 

1 . £6 

C . 25 

-C. 24 

0. 73 

Oo 03 

-0. 06 

-0.32 

0,60 

1 . 52 

0,33 

-C. 33 

0. 62 

0. 09 

0.0 

-C.9 1 

0.50 

1 . *3 

C. 1 C 

- 0 . 1 c 

0, 6 5 

0. 07 

-C . 07 

— 0.9 4 

0)50 

1 . 4A 

0.40 

-0. 4° 

c, se 

0.44 

0. 07 

-0.9 3 

0.50 

1 .40 

-0o C2 

-0o 08 

0. 89 

0. 79 

0.45 

-0.97 

0.50 

1 . 36 

-0.77 

0. 77 

0. 71 

0. 86 

0.70 

-1.00 

0.50 

1 . 32 

- 0 . s s 

0.93 

0, 57 

0. 72 

0. 09 

-1.00 

0,5 C 

1 o 28 

0 . c 

0. 0 

0, 0 

0. 0 

0. 0 

0.0 

o.eo 

1 . 56 

-0.17 

-0. 01 

0. 79 

0.0 

-0.25 

-0.65 

0.6 0 

1 a 52 

0.00 

-C.62 

0. 91 

0.46 

-0.20 

-0.62 

0.80 

1.48 

C. 19 

-0. 55 

C. 85 

0. 6 7 

Oo 13 

-0.8 7 

0,8 0 

1.44 

o. ce 

-0.56 

0, 87 

0.83 

0.38 

-0.85 

0.00 

1 .-.0 

-0. 19 

-0.54 

0o 91 

0. 90 

0.41 

-0.7 1 

o.eo 

1 , 36 

-0.66 

-O. 30 

C. 94 

0.93 

0 . 30 

-0.5 1 

0.80 

1.32 

-C. 6*4 

C. 34 

0, 94 

0.91 

-0.C8 

-0.78 

0.6 0 

1 .28 

- 0 o So 

0.6 9 

0. 90 

0. 70 

-0.53 

-0.9 2 

1.00 

1 . 56 

— 0.29 

-0. 1 1 

0. 61 

-0. 08 

-0.34 

-0.09 

1.C0 

1 o c 2 

-0.64 

-0. 30 

0. 66 

0. 3 3 

-0. 13 

-C.26 

l.CO 

1 a A a 

-0,46 

-0. 52 

0. 7 3 

0.63 

0.01 

-0.44 

1 .00 

1 , 44 

0. 06 

-C. 74 

0. 94 

0.82 

0.03 

-0.70 

1 .00 

1.40 

C. 29 

-0. 84 

0. 90 

0.90 

0.01 

-0.8 3 

I .00 

1 . 36 

- C . 06 

— C. 68 

0. 83 

0. 83 

0. 01 

-0.7C 

1.00 

1 . 32 

-0.50 

-C.43 

0, 77 

0.74 

-0. 07 

-C.56 

1.00 

1 .2e 

-0.69 

-0. 1 4 

0. 76 

0,62 

-0 . 24 

-0.59 

1.25 

1 . 56 

0.01 

-0.49 

Oo 4 9 

-0.4 0 

-0.52 

-0.4 3 

1.25 

1 . =2 

0.13 

-0. 56 

0, 5 0 

0. 07 

-0. 34 

—0 .64 

1.25 

1.4? 

0. 10 

- C. 8 5 

0, 75 

Oo 5 8 

-0. 32 

-C.5 1 

1 .2 = 

1 . 44 

0.44 

-C. 90 

0, 65 

0.82 

-0. 45 

-0. 76 

1.25 

1 . 40 

0.61 

- C. 8 8 

0, 6 6 

0. 85 

-0.44 

-0.90 

1.25 

1 . 36 

-C» 09 

-0. 57 

0. 52 

0.49 

-0. 13 

-C.76 

1.25 

1 . 32 

-C .56 

-0.42 

C, 37 

0.23 

-0. 24 

-0.4 9 

1 .25 

1.28 

— 0.74 

-0. 4 8 

0, 4 6 

-0. 08 

-0.54 

-0.0 8 

2.00 

1 . 56 

0.39 

- C. 66 

C. 4 2 

-0.4 l 

-0.42 

— C c 6 3 

2.00 

1 . 52 

0 o 59 

-o.ac 

0. 5 7 

-Oo 25 

-0.55 

-0.77 

2.00 

1 o 68 

0.52 

-0.81 

0. 5 5 

-Oo 26 

-Oo 54 

-0.7 1 

2.00 

1 O 44 

0.57 

— C. 84 

Oo 46 

0. 19 

-0. 30 

-0. 78 

2.00 

1.40 

0 o 4 3 

-0. 83 

0, 50 

0.37 

-Oo 29 

-0. 70 

2.00 

1 . 36 

-C . 54 

-Co 95 

-0. 03 

0. 06 

Oo 26 

0.28 

2.0 0 
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See text for explanation of tabular material. 



